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ABSTRACT 

Several lines of evidence have suggested that the the galaxy cluster CI 0024-1-17, an apparently 
relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. 
In this paper we present a high-resolution iV-body/hydrodynamics simulation of such a collision. 
We have created mock X-ray observations of our simulated system using MARX, a program that 
simulates the on-orbit performance of the Chandra X-ray Observatory. We analyze these simulated 
data to generate radial profiles of the surface brightness and temperature. At later times, t = 2.0 — 3.0 
Gyr after the collision, the simulated surface brightness profiles are better fit by a superposition of two 
/3-model profiles than a single profile, in agreement with the observations of CI 0024-1-17. In general, 
due to projection effects, much of the post-collision density and temperature structure of the clusters 
is not seen in the observations. In particular, the observed temperatures from spectral fitting are 
much lower than the temperature of the hottest gas. We determine from our fitted profiles that if 
the system is modeled as a single cluster, the hydrostatic mass estimate is a factor ~2-3 less than the 
actual mass, but if the system is modeled as two galaxy clusters in superposition, a hydrostatic mass 
estimation can be made which is accurate to within ~10%. We examine some implications of these 
results for galaxy cluster X-ray surveys. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (CI 0024+17) — methods: 
N-body simulations — X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies have proven to be interesting labo- 
ratories for studying the dynamics of the different kinds 
of matter in the universe. Early measurements of galaxy 
velocity dispersions in clusters demonstrated that most 
of the mass of galaxy clusters was in the form of a non- 
luminous component (jZwickylll937f ). Modern cosmologi- 
cal constraints (e.g. Turner 2000) suggest that this non- 
luminous material must be mostly non-baryonic. In ad- 
dition, X-ray observations of clusters of galaxies reveal 
that most of their baryonic mass is in the form of a hot, 
diffuse, X-ray emitting gas (the intracluster medium or 
ICM). The presence of these different types of matter 
(stellar, gaseous, and dark) in these systems provides not 
only insights into the nature of the clusters themselves 
but also the nature of the different components of matter 
and their dynamics. 

In the cold dark matter (CDM) picture of structure 
formation, galaxy clusters are the product of many merg- 
ers between galaxies, galaxy groups, and smaller clusters 
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of galaxies (|Davis et al.lll985f ). Observations of galaxy 
clusters indicate that this process is still ongoing in many 
systems. A famous example is the so-called "Bullet Clus- 
ter" (IE 0657-56), a merging system in which X-ray and 
weak gravitational lensing observations demonstrate a 
clear separation between a collisionless dark matter com- 
ponent and the intracluster medium (|Markevitch et all 
2002; Clowe et al. 2006). 

Another system that has recently attracted attention 
is CI 0024-1-17, a cluster at a redshift z = 0.395 ex- 
hibiting both weak and strong lensing. Early attempts 
at reconstru cting the mass profile of this system us- 
ing lensing (iTvson et all 119981 iBroadhurst et all 120001 : 
iComerford et al.ll2006D suggested a conflict with the pre- 
dictions of the standard CDM model due to the flatten- 
ing of the density profile in the inner regions of the clus- 
ter. Cosmological simulations assuming CDM indicate 
that galaxy cluster mass profiles should exhibit a nonzero 
logarithmic slope in the inner regions (e.g. Navarro, 
Frenk, & White 1997). This apparent discrepancy 
(and others) led to suggesti ons that the CDM paradigin 
would need to be m odified CSpergel fc Steinh ardtl l2000l : 
iHogan fc Dalcanton 2000; Moore et al. 2000), for exam- 
ple to include self-interaction of the dark matter. 
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ICzoske et all (|2001l I2002D demonstrated that the red- 
shift distribution of the cluster galaxies in CI 0024+17 is 
bimodal, with a large, primary component and a more 
diffuse, secondary component, separated from the pri- 
mary component in velocity space by w ~ 3000 km/s. 
They suggested that the system is composed of two clus- 
ters undergoing a high-velocity collision along the line 
of sight. They also performed a simulation demonstrat- 
ing that such a scenario reproduces not only the bi- 
modal redshift distribution but also the flattening of the 
central density pro file. X-ray obse rvations of the clus- 
ter using Chandra (fOta et al.l[200l and XMM-Newton 
(jZhang et al.ll2005l ) revealed that the surface brightness 
profile is better fit by a superposition of two ICM mod- 
els rather than one. They suggested on the basis of the 
isothermal temperature profile that the individual clus- 
ters of the system had returned to equilibrium after the 
collision and that consequently the encounter must have 
occurred several Gyr ago. 

Recentl y, a weak lensing analysis presented by 

iJee et all ()2007[ ) revealed a ringlike structure in the pro- 
jected matter distribution. They proposed that dark 
matter from the cores of the clusters had been disrupted 
and ejected from the systems by the collision, and they 
demonstrated that such features could be reproduced in 
a simulation of a collision of two pure dark matter halos. 
From this they suggested the current state of the system 
is ~l-2 Gyr past the pericentric passage of the cluster 
cores. We are addressing this issue in a separate paper 
(ZuHone et al. 2008, in preparation). iJee et all (|2007l ) 
also demonstrated that if the system is modeled as two 
ICM profiles in superposition, the mass determination 
based on hydrostatic equilibrium agrees with the result 
from their lensing analysis. 

If the merger scenario for CI 0024-1-17 is correct, it 
raises several important questions: how long after the 
collision would the clusters appear relaxed, and corre- 
spondingly when would a mass estimate based on hy- 
drostatic equilibrium yield an accurate measurement? 
What effect does viewing such a collision along the line 
of sight, with both cluster components in superposition, 
have on the observed density and temperature structure? 
In this paper we seek to provide answers to these ques- 
tions by simulating a similar high-speed collision between 
two clusters of galaxies. Most importantly, we include 
the dynamics of the cluster gas, which were not included 
in the previous simulations. 

Throughout this paper we assur ne a flat CDM co smol- 
ogy with h = 0.5, = 1.0, as in lOta et all (|2004f) . 

2. SIMULATIONS 
2.1. Method 

We performed our simulations using FLASH, a parallel 
hydrodynamics/ A^-body astrophysical simulation code 
developed at the Center for Astrophysical Ther monu- 
clear Flashes at the University of Chicago (Frvxei ret al.l 
[2003). FLASH uses adaptive mesh refinement (AMR), 
a technique that places higher resolution elements of the 
grid only where they are needed. In our case we are 
interested in capturing sharp ICM features like shocks 
and "cold fronts" accurately, as well as resolving the in- 
ner cores of the cluster dark matter halos. As such it 
is particularly important to be able to resolve the grid 



adequately in these regions. AMR allows us to do so 
without needing to have the whole grid at the same res- 
olution. FLASH solves the Euler equations of hydrody- 
namics using the Piecewise- Parabolic Method (PPM) of 
IColella fc WoodwardI (jl984) , which is ideally suited for 
capturing shocks. FLASH also includes an A-body mod- 
ule which uses the particle-mesh method to solve for the 
forces on gravitating particles. The gravitational poten- 
tial is co mputed usin g a multigrid solver included with 
FLASH (|Rickejl2008D . 

2.2. Initial Conditions 

The clu ster dark matter halos are init ialized using the 
NFW (.Navarro. Frenk. fc Whitelll997l) density profile , 
where we follow the method of iKazantzidis et al.l (|2006f l . 
The NFW functional form is used for r < r2oo- 



PDM(f) = 



r < 7-200 



(1) 



Outside of r2oo , to keep the mass of the halo finite and to 
avoid an unphysical distribution function, we implement 
an exponential cutoff that turns off the profile on a scale 
'"docayj ^ free parameter which we set to 0.1r2oo: 
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Here r2oo is the radius within which the mean mass den- 
sity is 200 times the cosmic mean, is the NFW scale 
radius, and C200 = ''200 /''s is the NFW concentration 
parameter. We require that at r2oo the profile and its 
logarithmic slope be continuous. This is achieved via the 
parameter k: 



(1 + 3C; 
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(3) 



Truncating the cluster NFW profiles at r2oo is unphys- 
ical. In real dark matter halos, the density profile is 
expected to remain NFW-like out to where the density 
of the halo reaches the average matter dens ity of the uni- 
verse (|Prada et al.ll2006l : iTavio et al.ll2008[ ). Single clus- 
ter tests we have performed (as well as those performed 
in other studies, e.g., Ricker and Sarazin 2001) demon- 
strate that the truncated density profile is stable for the 
length of the run, with minor steepening of the profile 
near the cutoff radius. 

The truncation also reduces the halo mass by approx- 
imately 25%. Assuming the clusters fell in from infinity, 
the larger halo masses would accelerate the halos to an 
impact velocity approximately 12% larger than the value 
we use. After the collision, the corresponding decelera- 
tion due to gravity would also be higher. The effect of the 
larger halo masses on the deceleration due to dynamical 
friction (/dyn on pA-P/v"^) is more comphcated to deter- 
mine, since although the density in the outskirts of the 
halos is higher (thus increasing the dynamical friction), 
the velocity of clusters would also be higher due to the 
increased acceleration (thus decreasing it). The higher 
infall velocity of the clusters and the increased mass out- 
side r2oo might also lead to some increase in the dissipa- 
tion of the shock energy, and therefore in the heating of 
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TABLE 1 
Initial Cluster Parameters 



Cluster 


M200 (Mq) 


»"200 (kpc) 


C200 


ric (cm ^) 


Tc (kpc) 


1 


6.0 X IQi* 


1739.79 


5.0 


0.020 


198.83 


2 


3.0 X IQi-* 


1380.87 


7.0 


0.062 


98.65 



the gas in the outskirts. However, in general any such ef- 
fects are inaccessible to current X-ray observations, due 
to the low densities in the cluster outskirts. We conclude 
that the effects due to truncation of the cluster density 
profile are not likely to be significant for the purposes of 
this study. 

To initialize the gas density we choose a Burkert profile: 



Pc 



{l + (r/re)2)(l + r/re) ' 

which orig inally (|Burkertl Il995l) was given as a fit- 
ting function for dark matter profiles of dwarf galax- 
ies but is also a good fit to the gas density profiles 
of non-cooling-core clusters (A. Kravtsov 2007, pri- 
yate communication). The mor e traditional /3-model 
(|Cavaliere &: Fusco-Femianolll976( ) is a poorer fit to the 
gas density profile at large r in both simulations and 
observations (IVikhlinin et al.ll2006l : iHallman et"all2007l : 
iNagai et al.ll2007p (However . we follow lOta et al.l ()2004f ) 
and IZhang et al.l ([20 05) and fit /3-model profiles to the 
X-ray surface brightness profiles in our simulated obser- 
vations.) 

The gas density is also fitted to an exponential taper at 
r2oo, extending to the radius at which it equals the mean 
baryonic density of the universe. For the gas profiles, we 
choose the core radius to be roughly half the scale ra- 
dius of the DM profile (|Ricker fc Sarazinll2001[ ). We fix 
the gas mass fraction at r2oo to 0.12 in line with observa- 
tions of real clusters (e.g., Vikhlinin et al. 200 6, Sander- 
son et al. 2003), and the measurement of lOta et al.l 
(l200l . 

From these parameters the central gas densities pc can 
be determined. The equation of state for the gas is an 
ideal-gas equation of state with 7 = 5/3 and mean atomic 
mass ^ — 0.592. We determine the pressure, tempera- 
ture, and internal energy profiles by assuming our func- 
tional forms for dark matter and gas density and numeri- 
cally integrating the equation of hydrostatic equilibrium. 
The values used for the cluster parameters are given in 
Table 1. 

FoUowing lCzoske et al.l ()2002[ ) and lJee et all (|2007f) we 
assume a mass ratio of 2:1 for the clusters. In our sim- 
ulation the clusters are initialized so that their centers 
are separated by the sum of their respective r2oo val- 
ues (approximately 3 Mpc), and they are given an ini- 
tial relative velocity Wrei = 3000 km/s. This value is 
chosen because it is the inferred relative velocity of the 
two redshift compo nents of CI 0024-1-17 seen in the red- 
shift histogra ms of iCzoske et"all (l200l. T he results of 
iCzoske et ahl ()2002[ ) and iJee et alTirsOOTD suggest the 
clusters have already com pleted their first pericentric 
passage, and ICzoske et al.l |2002) show that in their col- 
lision simulation an initial velocity of ^3000 km s~^ is 
needed to reproduce the observed distribution of red- 
shifts after the collision. This value for the velocity is 



approximately 2w^/ for the clusters, where Vff is the free- 
fall velocity from infinity. Such an initially high relative 
velocity might be difficult to achie ve in a ACDM universe 
(see, e.g. iHavashi fc White! 1200^ . However, our choice 
is motivated by our desire to match the relative inferred 
velocity of the cluster components of CI 0024-1-17. 

We refine the adaptive mesh using the second deriva- 
tive of density, temperature, and pressure to capture the 
sharp features in the gas dynamics, and using the mat- 
ter density to resolve the cores of the clusters. For our 
box size of 14.29 Mpc we achieve a finest resolution of 
Ax = 13.96 kpc. 

2.3. Simulated Observations 

To make a meaningful comparison with the X-ray ob- 
servations of CI 0024-1-17, we construct mock X-ray ob- 
servations of the simulation output. To do this we use 
MARX, a suite of programs created by the MIT/CXC 
group that simulates the on-orbit performance of the 
Chandra X-ray Observatory MARX provides detailed 
ray-tracing simulations of astrophysical sources as ob- 
served by Chandra. In this paper we use MARX to gen- 
erate simulated surface brightness profiles and spectra. 

For input into MARX we create projected flux maps 
from the FLASH sim ulation data. For each FLASH zo ne 
the MEKAL model l|Mewe. Kaastra. fc Liedahllll995D is 
used to assign a spectral emissivity. For these simula- 
tions zero metallicity has been assumed. The spectra are 
projected along the appropriate line of sight, and fiux 
maps are generated for photons in 198 separate energy 
bins from 0.1-10.0 keV, giving an energy resolution of 
AE = 50 eV, compared to the AE « 100 eV resolution 
of Chandra. 

For our simulated Chandra observations we use the 
simulated ACIS-S detector and HRMA mirror system, 
as in the real Chandra observation of CI 0024-1-17. Each 
observation is ^40 ks of simulated exposure time, compa- 
rable to the Chandra observation of CI 0024-f 17 (39 ks). 
The flux maps serve as a distribution function for the 
simulated photons and determine the direction cosines, 
energies, and timesteps for the photons in MARX. For 
simplicity, effects such as foreground galactic absorption 
and point source contamination are not included. How- 
ever, we do include a diffuse X-ray background compo- 
nent. Output is given as a FITS file that can be read 
and analyzed in the same way as real Chandra observa- 
tions. To analyze our simulated data, we use CIAO 3.4 
and XSPEC 12.3. Details of our procedure for generat- 
ing these observations and a verification study are given 
in the Appendix. 

3. RESULTS 

3.1. Qualitative Description 

Figures 1 and 2 show slices of density and temperature, 
respectively, through the z — coordinate plane at dif- 
ferent times in the simulation (the x axis is the collision 
axis). From t = 0.0-0.9 Gyr the clusters approach each 
other and a shock front forms. At t — 1.0 Gyr the core 
of the smaller cluster collides with the core of the larger 
cluster, driving a shock wave forward into the larger core 
and displacing the larger core's gas from the collision axis 

^ http://space.mit.edu/ASC/MARX/ 
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Fig. 2. — Slices of gas temperature through the z ■ 
scale is logarithmic. Each panel is 10 Mpc on a side. 



coordinate plane for times t = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 Gyr. The color 



into streams of cold gas. There is also a stream of dense, 
cold gas that is pulled directly between the two dark mat- 
ter halos. Later on, the heads of these streams of gas, 
along with hotter gas from between the clusters, begin 
to fall into the potential of the larger halo and by t — 2.0 
Gyr have fallen in completely. As for the gas core of the 
smaller cluster, it is penetrated shortly after the collision 
by a reverse shock, but this shock weakens as it traverses 
the core and so the latter remains cool, in agreement with 
similar investiga tions of high- velocity clu s ter m ergers by 
iTakizawal (|2005[) and lMilosavlievic et all (|2007D . In this 
process some of the gas of the smaller cluster is also 
ram-pressure stripped, and a contact discontinuity forms 
that is similar to the one observed in the 1E0657-56, and 



in the simulations of the latte r (jSpringel fc Farrai1l2007l : 
iMastropietro &: Burker"tll2008[ ). 

After penetrating the larger's gas core the smaller core 
is forced to lag behind its dark matter core by ram pres- 
sure until it encounters regions of lower density {t > 1.2 
Gyr). At this point the gas falls back toward the dark 
matter core but overshoots the potential and begins to 
slosh back and forth inside it. At late times the gas of 
the clusters is still settling into their respective dark mat- 
ter potential wells, and the temperature structure of the 
gas is still complicated. Due to the high initial veloc- 
ity of the collision, the clusters are in an unbound orbit 
and by the end of the simulation the separation between 
their respective dark matter halo centers is still increas- 
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Fig. 3. — Mock 40 ksec raw counts image of the simulation viewed along the line of centers of the clusters in the 0.5-5.0 keV band. The 
image has been smoothed using a Gaussian with a fixed kernel radius of 3 pixels. Color scale units are counts/pixel. The time is t = 3.0 
Gyr, ~2 Gyr after the collision. 




Fig. 4. — Plots of surface brightness vs. 
times t = 2.0, 2.5, and 3.0 Gyr. 



projected radius for 



ing. This is in contrast to many previous works involving 
mergers (e.g., Ricker & Sarazin 2001; Poole c t al..,2006.) . 
which assumed an initial velocity v w Vff. 

3.2. Mock X-ray Observations 

For the X-ray observations we view the system along 
the line of centers. Since we do not know the current 
state of the actual collision, we set each individual obser- 



vation as if it were being observed at a redshift z = 0.395. 
At this redshift for the assumed cosmology 1" = 6.39 kpc. 
Figure 3 shows an example of a smoothed mock X-ray 
image. For each observed time we extract a radial surface 
brightness profile and a radial temperature profile. 

Cluster surface brightness profiles are extracted for the 
energy range 0.5-5.0 keV, within a circle centered on the 
surface brightness peak with a radius of 400" (~ 2500 
kpc). Figure 4 shows the surface brightness profiles for 
times t = 2.0, 2.5, and 3.0 Gyr. The profiles are fitted 
with the /3- model (Cavaliere & Fusco-Femiano 1976): 



Sx{r) = So 



1 



-3/3+1/2 



(5) 



We fit with a single /3-model and a sum of two /3-models. 
The Y^-st a tistic is us ed to fit t he profiles. We follow 
lota et all (|2004D and iJee et all (2007) in fixing one of 
the P parameters to unity, as there exists a strong de- 
generacy between this parameter and its corresponding 
core radius r^- We include the background constant as a 
free parameter in the single /3- model fit, but freeze this 
value for the double /3-model fit. The fitted profiles for 
the t = 3.0 Gyr case are shown in Figure 5. The values 
of the fitted parameters for epochs t = 2.0, 2.5, and 3.0 
Gyr are given in Tables 2 and 3. Errors are quoted at 
the 90% confidence level. 

The projected temperature profiles for t — 2.0, 2.5, 
and 3.0 Gyr are shown in Figure 6. The profiles were 
generated by extracting spectra from 4 annular regions in 
the 0.5-7.0 keV band centered on the peak of the surface 
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t = 3.0 Gyr, Double /S-model fit, /Sj frozen 




I = 2.0 Gyr h 
t = 2.5 Gyr ^ 
t = 3.(VGyr ^ 



100 200 300 400 500 600 700 
r (kpc) 



Fig. 6. — Plot of temperature profile vs. projected radius for 
times t = 2.0, 2.5, and 3.0 Gyr. 



brightness profile. Each is fitted with a MEKAL model, 
using the x'^-statistic with grouping to ensure at least 15 
photons per bin. The radial extent of the temperature 
profiles is 125" (« 800 kpc), which is slightly greater 
than the radial r ange for w h ich th e temperature profile 
was extracted in lOta et all (|2004[ ) . Single-temperature 
fits to the spectrum of the whole cluster in the 0.5-7.0 
keV band for t = 2.0, 2.5, and 3.0 Gyr are given in Table 
4. 

4. DISCUSSION 

4.1. The Surface Brightness and Temperature Profiles 

The feature in the rad ial p rofile of C I 0024+ 17 pointed 
out by lota et all (Hooi) and lJee et all (|2007l ) was that a 
double /9-model fit is a better fit to the data, as a single /3- 
model fit underestimates the central surface brightness. 
The surface brightness profile of CI 0024-1-17 and the sin- 
gle and double /3-model fits to the profile are shown in 
Figure 7. The difference between the fits is significant; 
there is a difference in the x^-statistic of ~60 for a dif- 
ference of two degrees of freedom. 





Data I--+-H 




Single p-model fit 




Double p-model fit 









:* 



















Fig. 7. — CI 0024+17 X-ray surface brightness profile and sin- 
gle and double /3-model fits. Reproduced with permission from 
iota et al.l ll200l) . 



In our mock cluster observations, we find that at later 
times, it is also possible to make a distinction between 
the two surface brightness components of the two clusters 
from the difference in the fitted x^-statistic (see Tables 
2 and 3). At the t = 2.0 Gyr epoch, Ax^ «45 for a dif- 
ference of only 2 degrees of freedom, and at later times 
is even larger (Ax^ Ri200-260 for a difference of two de- 
grees of freedom). We also note that we can use the 
Ax^-test to distinguish between the two clusters in our 
initial conditions, which we also show in Tables 2 and 3. 

In CI 0024-1-17, the fitted spectral temperature of 
T = 4.47 keV seems low for a system undergoing a col- 
lision or merger. In our mock observations we note a 
similar phenomenon; the fitted temperatures in the mock 
observations are significantly lower than what might be 
expected given the temperatures observed in the simu- 
lation. In our simulated system, the hottest cluster gas 
is in the larger cluster and is at a temperature of ~ 5-6 
keV at later times. However, the measured temperatures 
in the mock X-ray observations barely reach ^ 3 keV. A 
close look at the temperature distribution in the Simula- 
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TABLE 2 

Fitted Cluster Parameters for a Single /3-model 



t (Gyr) 


Tc (kpc) 


So (cts s -"^cm ^arcsec ^) 


/3 


xVd.o.f. 


0.0 
2.0 
2.5 
3.0 


172.04l'[;^| 

9*^^ 91 +2.24 


(3.8lt0;«2) X 10- 
(1.42l0;°|) X 10- 
(1.25l° °|) X 10- 
(1.43l0;0|) X 10- 


-6 
-7 
-7 
-7 


n 7QO+0.001 
n 7i4+°-°i'' 

n Qfi4+0 009 
U.MD^_0.03l 

n yoQ+O.OOS 

U. / JO_Q 027 


407.61/195 
305.14/195 
491.04/195 
423.58/195 



TABLE 3 

Fitted Cluster Parameters for a Double /3-model 



t (Gyr) rci (kpc) 5o,i (cts s -"-cm 



Pi 



(kpc) 5o,2 (cts s 



xVd.o.L 



0.0 
2.0 
2.5 
3.0 



+3.16 
0.32 
+34.41 
-17.45 
+9.04 
5.36 

82.12t«;i2 



89.24 
39.91 
79.28 



{3.73l°'0^) X 10- 

+0.17-, ^ 



(8.18 
(2.26 
(2.51 



-0.17/ 

+-!«) X 10- 



1.0 (F) 
1.0 (F) 
1.0 (F) 
1.0 (F) 



204.65 
270.00 
282.00 
316.24 



1.32 
4.00 
+6.34 
4.78 
+ 11.97 

5.57 
+26.16 
4.89 



(7.03 
(Lll 
(7.70 
(4.94 



tlfr) X 10- 


7 


0.981 


t^ol) X 10- 


7 


0.970 


+_rA) X 10- 


-8 


0.970 


tllt) X 10- 


-8 


0.980 



-0.014 
-0.010 
-0.025 
-0.021 
-0.021 
-0.015 
-0.049 



260.01/193 
245.83/193 
226.75/193 
205.86/193 



Fig. 
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33.00 



- 9.20a 



2.569 



0.7168 
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Slice through the 2 = coordinate plane of the simulation of temperature at t = 3.0 Gyr, ~2 Gyr after the collision. The 



colorscale is logarithmic. Density contours are logarithmically spaced by 2.5x, between n : 
Mpc by 6 Mpc. 



10- 



■ 10 cm . The size of the figure is 14 



TABLE 4 
Fitted Average Spectral 
Temperatures 



t (Gyr) 


Tspcc (keV) 


X^/d.o.f. 


2.0 
2.5 
3.0 


^•"'-0.16 

QQ~l~t).23 

^•"■^-0.18 

r, 71+0.20 

^•'J^-0.19 


140.62/177 
147.98/174 
128.07/167 



tion itself (Figure 8) reveals that the highest-temperature 
gas is confined to the central region of the larger clus- 
ter, and the denser gas of the smaller cluster is colder, 
around 3 keV. In projection along the line of sight, the 
lower-temperature, denser gas "washes out" the higher- 
temperature, less dense gas, and the complicated, two- 
component temperature structure of the clusters is lost 
because of the strong density dependence of the X-ray 



200 300 400 500 
r (kpc) 



600 700 



Fig. 9. — Radial profile of mass- weighted temperature at t = 3.0 
Gyr, ~2 Gyr after the collision. 



emission (oc p^). Figure 9 shows the profile of the mass- 
weighted temperature over the same range as Figure 6 for 
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TABLE 5 

Estimated and Exact Masses at R = 35" 



t (Gyr) Mx,I3(Mq) Mx,2/3(M0) A^ctualCMe) Mx./s/A^actual Mx,2f(/Mactual 

2.0 (e.oalg'ae) X 10" (1.7llE! o^) x lO" 1.74 x lO" 0.35 0.98 

2.5 (7.50l°;|^) X 10" i'^-'^'^tom) ^ 1-65 x 10" 0-45 1.05 

3.0 (6.57lg;^;^) X 10" (1.62j;f,:J°) x 10^* 1.55 x 10^* 0.42 1.05 



3 I 1 ' ' ' ' 1 

100 200 300 400 500 600 

r (hgo"' kpc) 

Fig. 10. — CI 0024+17 Temperature profile. Average profile is 
shown as a dashed line at T = 4.47 keV. Reproduced with permis- 
sion from iota eFaP ll200i ). 

t = 3.0 Gyr. The lack of temperatures higher than 3 
keV demonstrates that the denser, colder gas provides 
the dominant contribution to the X-ray emission. 

The single-temperature fits are in general agreement 
with our projected temperature profiles. Attempting 
to fit the cluster temperature using a two-temperature 
model resulted in no significant decrease in the x^- 
statistic from the single-temperature fits, so on the basis 
of this analysis we cannot distinguish the two different 
average temperatures of the two clusters. 

At late times there is still some moderate variation in 
the projected temperature profile with time, indicative 
of the continuing evolution of the system, though the 
average temp erature of the cl uster is roughly constant 
with time. In lOta et al.l (120041 ) . the observed isothermal 
temperature profile of CI 0024-1- f 7 (reproduced in Figure 
10) was taken as evidence that the merger was not recent 
(i.e., it was more than a few Gyr ago). Our results show 
that if viewed along the line of sight ~2 Gyr after the 
merger, the temperature structure of the system can ap- 
pear more regular in projection than it actual ly is (e.g., 
the structure seen in Figure 2). If the results of ljee et alj 
(jlOO?) regarding a ringlike dark matter structure indi- 
cate the clusters are being viewed along the line of sight 
shortly after the merger, our results show that the sys- 
tem can still appear relatively relaxed. This agrees with 
previous results of mock X-ray observations of simula- 
tions of galaxy cluster mergers (e.g., Ricker and Sarazin 
2001, Poole et al. 2006). 

Finally, our observed spectral temperature from our 
simulations ofT ~ 2.6 — 2.8 keV is somewhat lower 
than the T ~ 4 .47 keV for CI 0024-^17 measured by 
Ota et al I (I2004D or the T - 4.25 keV measured by 
Jee et al.l (j2007). If the collision scenario for CI 0024+17 
is correct, this indicates that the initial conditions for 
our simulation are not identical to the pre-merger con- 



ditions for the real cluster. A higher impact velocity or 
a change in the masses of the clusters could account for 
the temperature difference. 

4.2. The Reliability of Hydrostatic Mass Estimates 

Though the temperature and surface brightness pro- 
files are undergoing considerable evolution even at late 
times, it is still relevant to ask how accurate a mass es- 
timate would be under the assumpti on of hyd r ostatic 
equi librium. As previously noted, lOta et ahl ()2004D 
and iJee et all (2007) made estimates of the mass of 
CI 0024-1-17 using single and double /3- model fits to the 
surface brightness profile together with an isothermal 
temperature profile. Although our temperature profile 
is not strictly isothermal, existing methods for tempera- 
ture deprojection assume spherical symmetry, so we do 
not attempt such a deprojection here. If we assume a 
cluster temperature Tgpec given by the spectral temper- 
ature fit to the entire cluster, we can arrive at a rough 
estimate for the mass of the system under the assumption 
of hydrostatic equilibrium. 

For a /3-model fit to the cluster gas, the projected clus- 
ter mass within a cylindrical volume for a given tem- 

J ierature, B, ari d core radius Tc can be estimated by 
Ota et al.ifl998l: iJee et al.l[200l 

Mxp{R) - 1.78xl0"/3 ( J-) ( ^) ^/''^ Mq , 

(6) 

where R is the projected radius, T is the gas tempera- 
ture and (3 is the /3-model index. We compute the esti- 
mated mass within the arc radius of the lensed galaxy 
in CI 0024-hl7, rare = 35'7223.65 kpc. For the double 
/3-model fits we assume the same temperature for both 
cluster components, and we add the mass contributions 
from the two fits together. Table 5 shows the estimated 
masses from the single and double /3-model fits for times 
t = 2.0, 2.5, and 3.0 Gyr. Errors are quoted at the 90% 
confidence level. It is clear from the table that assuming 
a single-cluster model for the system will underestimate 
the projected mass by a factor ~2-3, whereas assuming 
a double-cluster model will estimate the projected mass 
to within ~10%. This is in agreement with the results of 
IJee et al.l (l2007l) . 

4.3. The Effect of A Slightly Off- Axis Collision 

In our simulated head-on collision, the gas in the clus- 
ters is significantly disrupted and by the end of the sim- 
ulation is still settling into the potential wells created by 
the dark matter. If instead of colliding head-on the en- 
counter were allowed to be slightly off-axis, this would 
result in less disruption of the cluster gas cores. Presum- 
ably the gas would relax more quickly and an even more 
accurate mass determination could be made under the 
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assumption of hydrostatic equilibrium. In addition, in 
a slightly ofF-axis collision viewed along the line of sight 
the two cluster components would be more clearly dis- 
tinguished in the X-ray emission. Two-component fits to 
both the X-ray surface brightness and the temperature of 
the different components would enable a more accurate 
mass determination. This is a large parameter space of 
simulations that will be explored in future investigations, 
but is beyond the scope of this paper. 

4.4. Implications for X-ray Surveys 

These results have important implications for X-ray 
surveys. Neighboring galaxy clusters viewed in superpo- 
sition, whether undergoing mergers or not, may bias tem- 
perature measurements based on spectral fitting. This 
can have an effect particularly on X-ray-based scaling re- 
lations. For example, previous simulated observations of 
merging clusters have demonstrated that for mergers oc- 
curring along the line of sight, the clust ers appear to have 
higher temperatures fo r their mass (jRicker fc Saraziiil 
I2OOII : iPoole et ani2007D . Our investigation shows that 
clusters viewed in superposition can make a cluster ap- 
pear colder than it actually is if the colder cluster's gas is 
significantly denser. In addition, discrepancies between 
hydrostatic and lensing mass estimates may in some cases 
be attributed to the existence of multiple components 
viewed along the line of sight, as is apparent in the case 
of CI 0024-1-17. Detailed analysis of X-ray observations 
would be required to discern such superpositions. 

One possible way to correct for and identify such su- 
perposition effects is by using the optical redshifts of 
the cluste r galaxies to i dentify sepa rate cluster compo- 
nents, as iCzoske et all (I2OOII [200I) did in the case of 
CI 0024-1-17. To reliably identify such separate clus- 
ter components requires accurate redshifts. Because 
spectroscopic redshifts are only feasible for the nearest 
or brightest cluster galaxies, cluster surveys must rely 
on photometric redshift e stimators. Large optical sur- 
yevs such a s zCOSMOS (|Lillv et al.l 120071 ). Combo-17 
(IWolf et al.ll2003D. and the upcoming Dark Energy Sur- 



vey ( The Dark Energv Survev Collal3orationll2005( ) have 
photometric redshift errors cr^ ~ (0.02 — 0.05) (1 -I- z). 

If the only redshift determination for CI 0024-1-17 were 
via photometry, would the separation between the two 
components still be discernible? To answer this question 
we constructed a mock "galaxy" catalog from our sim- 
ulation. We chose dark matter particles from the two 
clusters as proxies for the cluster galaxies (120 from the 
smaller cluster, 240 from the larger) and assumed a value 
for CTz- We constructed redshift histograms for the re- 
sulting galaxy distributions and fitted them to both a 
Gaussian distribution and a sum of two Gaussian distri- 
butions: 

fs^ngle = AoC' f'"^" ^^'^0 (7) 



(8) 



The fitted parameters for these distributions for varying 
(Tz , given a choice of mock galaxies from our simulation, 
are shown in Table 6. For low values of CTz, the mean 
values of the two redshift distributions are statistically 
distinguishable, but for higher values they are not. The 
histograms shown in Figures 11 and 12 demonstrate that 
for az ~ 0.005 — 0.01 the skewness of the distribution of 



galaxies is evident, but when the redshift error increases 
to <7z ^0.02, the resulting galaxy distribution is indistin- 
guishable from that of a single cluster. 

At a redshift of z « 0.4 (CI 0024-1-17), the expected 
photometric redshift errors from optical surveys would 
be too large at this time to distinguish between the two 
cluster components (cz ~ 0.03 — 0.07), and at a redshift 
of z sa 1 it would be worse (cz ~ 0.04 — 0.1). Identifica- 
tion of line-of-sight mergers and collisions from redshift 
distributions likely will be restricted to clusters for which 
we have spectroscopic redshifts available. 

Our X-ray fitting results might be taken to suggest 
that X-ray surface brightness profiles could be used in 
place of galaxy redshift distributions to detect clusters 
seen in superposition. To test this hypothesis, we placed 
our simulated clusters at the epoch t — 3.0 Gyr (^2 Gyr 
after the collision, where we can clearly distinguish both 
cluster components) at two redshifts, the current red- 
shift of z = 0.395 and a redshift of z = 1.0. The same 
analysis of the surface brightness profiles as above was 
performed for exposure times of 10 ks, 40 ks, and 240 
ks, the former being a more typical exposure time for 
a given cluster in an X-ray survey. We focus here on 
the difference in the between the single and double 
/3-model fits. The resulting Ax^ for the fitted surface 
brightness models corresponding to the different expo- 
sure times for the different redshifts are shown in Table 
7. At z = 0.395, it is possible to use the Ax^ between 
the two models to distinguish between the two clusters 
for all exposure times. At z = 1.0, for an exposure time 
of 10 ks, the Ax^ is ~ 14 for a change of 2 degrees of 
freedom, and for 40 ks Ax^ ~ 70. Therefore, we find that 
at exposure times more typical of X-ray surveys, even at 
high redshift we can distinguish between the two cluster 
components, though with more difficulty. This is poten- 
tially a promising avenue for identifying merging clusters 
and clusters in projection. However, the ability to iden- 
tify merging clusters requires that the two clusters have 
surface brightness profiles that are distinct enough that 
the likelihood-ratio test can distinguish between them. 
For example, earlier on in the simulation, at the epoch 
t — 2.0 Gyr (w 1.0 Gyr after the collision, where we 
cannot distinguish between the two clusters as clearly), 
we find that the test can distinguish between the two 
components for our fiducial exposure time and redshift, 
but not at lower exposures and higher redshifts (Table 
7). Thus, determining exactly the physical circumstances 
under which the likelihood ratio test is the most powerful 
way to identify merging clusters and clusters in projec- 
tion requires further study. 

5. CONCLUSIONS 

The galaxy cluster CI 0024-1-17 is thought to be a col- 
liding system of galaxy clusters that is viewed along the 
collision axis. We have performed a simulation of a high- 
velocity head-on collision of galaxy clusters and created 
mock X-ray observations, viewing the clusters from the 
same direction. Previous investigations of this system 
had focused on the coUisionless dynamics of the galaxies 
and the dark matter; we also include gas in our simu- 
lation. The mock X-ray observations of our simulation 
indicate that the X-ray emitting gas is still undergoing 
moderate evolution at times '^1-2 Gyr after the collision. 
However, much of the complicated structure in density 
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TABLE 6 

Fitted parameters for Redshift Histograms with Varying ctz 



0-2 


HO 


TO 




Ml 


cri 


M2 


U2 


xVd.0.f.d„„(,(e 


0.005 
0.01 
0.02 
0.08 


0.9975 
0.9958 
0.9936 
0.9521 


0.0040 
0.0086 
0.0175 
0.0634 


12.6/10 
31.68/9 
8.64/9 
6.65/7 


1.002 
0.9998 
0.9947 
0.9733 


0.0040 
0.0086 
0.0175 
0.0634 


0.9959 
0.9889 
0.9742 
0.9260 


0.0052 
0.0089 
0.0126 
0.0794 


5.88/7 
14.7/6 
9.84/6 
3.56/4 



0.14 
0.12 



0.06 
0.04 
0.02 



1.005 1.01 



ol 6 

0.96 0.97 



Fig. 11. — Number of redshifts per bin vs. redshift, az = 0.005 and (Jz = 0.01. 



1,02 1.04 1.06 



0.6 0.65 0.9 0.95 1 1.05 1.1 1.15 1.2 



Fig. 12. — Number of redshifts per bin vs. redshift, az = 0.02 and az = 0.08. 



and temperature is obscured due to projection effects. 
The X-ray surface brightness profile is better fit by a su- 
perposition of /3-models than a single /3-model at times 
t = 1-2 Gyr after the collision, corresponding to the ob- 
served situation in CI 0024-1-17. The cluster gas in the 
center appears colder than it really is due to the projec- 
tion of denser, colder gas along the line of sight of the 
hotter gas. The temperature profile exhibits only mod- 
erate evolution with time, though the temperature dis- 
tribution of the clusters in the simulation has significant 
structure. Though a mass estimate of the cluster as- 
suming one cluster component in hydrostatic equilbrium 
underestimates the mass by a factor of ~2-3, a mass es- 
timate based on the assumption of two galaxy clusters in 
superposition comes close (~10%) to the actual projected 
mass of the system within the arc radius Taro in agree- 
ment with the mass measurements made of the cluster 
CI 0024-1-17 under the same assumptions. These results 
may be valuable when looking at clusters discovered in 



X-ray surveys. With current photometric redshift errors 
it will not likely be possible to distinguish such line-of- 
sight collisions from single clusters discovered in X-ray 
surveys at high redshift by identifying separate galaxy 
concentrations. However, our results show that it may 
be possible to distinguish merging clusters and clusters in 
projection from single clusters by testing multiple model 
fits against single model fits, though this will become 
more difficult at higher redshifts. 
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TABLE 7 

Ax^ FOR Different Exposure Times and Redshifts 







t = 3.0 Gyr 






t = 2.0 Gyr 






z = 


0.395 z = 


1.0 


2 = 


0.395 


z = 


1.0 


10 
40 
240 


177.69/96 
425.31/195 
1497.05/195 


- 113.21/94 99.20/96 - 

- 203.44/193 242.51/195 - 

- 439.33/193 770.37/195 - 


86.05/94 

173.29/193 

287.46/193 


96.41/96 
305.14/195 
428.73/195 


- 94.40/94 

- 245.83/193 

- 245.83/193 


115.42/96 - 
217.34/195 - 
280.26/195 - 


111.74/94 

209.53/193 

228.22/193 
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APPENDIX 

MOCK OBSERVATION GENERATION AND VERIFICATION STUDY 

Our mock X-ray observation generation procedure consists of two steps; we first use the FLASH input of density 
and temperature to c reate a flux m a p in s ky coordinates for a range of energies. In doing this we foUow closely the 
procedure outlined in lGardini et alj (|2004[ ). Secondly, we use this flux map as a "user source" in MARX to generate 
the photon energies, positions, and times for our simulated observation. 

The FLASH AMR grid is regridded to a uniform grid at the highest resolution in the simulation. We then choose 
a direction along which the cluster is observed and the physical quantities are projected. For each cell in our dataset 
the emissivity of the plasma is given by 

e = neni/A(r,Z) (Al) 

where rig and nn are the electron and hydrogen densities, respectively, and A(r, Z) is the power coefficient which 
depends on temperature and metallicity, which are assumed constant over one FLASH cell size. This coefficient is 
calculated using the MEKAL model. Using this relation, the photon luminosity (photons cm~^ keV^^) at a given 
energy hv is given by 

L2= [ eldV' = KZEM (A2) 
Jv 

where the quantity EM = jy n^nndV' is the emission measure. The measured flux of photons at a specific energy at 
a redshift z is given by 

where D/, is the luminosity distance for the given redshift and cosmology. 

This flux map is then used to generate photons in MARX using the "user source" implementation. The map is 
taken as an input and used as a distribution function for the photons to determine their position in sky coordinates, 
energy, and time of detection. In the cases considered here, the FLASH zone size is larger (typically by a factor of 
~4) than the size corresponding to a Chandra pixel. Within this FLASH zone size the photon positions are uniformly 
distributed. After the events file is generated, we re-bin the image by this factor to perform our spatial analysis. 

For the spatial analysis, we construct an exposure map for the ACIS-S chips with the input spectrum of the cluster 
as a set of weights, and when performing the fit we convolve the model with a 1-D PSF. For the spectral analysis, we 
construct weighted RMFs and ARFs (Redistribution Matrix Files and Auxiliary Response Files), and use them to fit 
a background-subtracted spectrum. Our spectra are grouped so that each bin has at least 15 counts. 

In order to ensure the accuracy of our MARX simulations, we have performed a few simple verification tests. We 
have constructed two models: one, a gas sphere with a constant temperature and a /3-model density profile, and the 
other, two gas spheres, each with different temperatures and /3-model parameters. We then subjected them to the 
same analysis as our simulated clusters in FLASH. In the single-cluster case, we verify that the surface brightness fits 
we recover the core radius Tc, /3-parameter, and central surface brightness Sq. Also, we verify in the spectral fit that we 
recover the model temperature for the fit to the whole cluster, the normalization of the spectrum, and the isothermal 
temperature profile. For the two-cluster case, we verify that the two cluster components can be distinguished and that 
the corresponding /3-model parameters are recovered (as in the simulated clusters, we find that we cannot distinguish a 
single spectral temperature model from a sum of spectral temperature models). For these tests we consider the model 
verified if we can recover the model parameters within the Icr errors. 

For the spatial fitting, we find that as we increase the exposure time of the observation, systematic effects become 
more important. At an exposure time of texp = 240 ks we recover the fitted parameters for both the single and double 
/3-model cases. This corresponds to the highest exposure time for our simulated clusters. These results are shown in 
Tables 8 and 9. However, as we go to higher exposure times, we find that we cannot recover the input models. Table 
11 shows the number of parameters which are not recovered and the corresponding reduced x^-statistic for increasing 
exposure time. For an exposure time of 1 Ms, we can verify the single /3-model but cannot verify the double /3-model, 
and for an exposure time of 4 Ms we cannot verify either. 

Possible systematic effects may include the modeling of the PSF (which is energy and spatially dependent), ACIS 
read-out errors and the precise dependence of the exposure map on the input spectrum. We note that for the fits 
that we cannot recover the parameters the resulting x^-statistic is large, indicating that the systematic errors are 
dominating the statistical errors. We also note that even for these cases we find the differences in the fitted parameters 
and the true parameters are at most on the order of a few percent. 

For the spectral fitting we find that up to a high exposure time of 4 Ms that we recover the input temperature and 
the normalization of the spectrum, as well as the isothermal temperature profile. The results of the fit are shown in 
Table 10 and the temperature profile is shown in Figure 13. 
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TABLE 8 

Fit to Single Cluster Test, te^p = 240 ks 



Type 


r-c (kpc) 


P 


So (cts s ^cm ■^arcsec ^) 




Model 
Fitted 


100 


1 

1 010+°°°^ 

1.U1U_(, Q21 


8.26 X 10-9 
(8.18^°;^^) X 10-9 


197.92/195 



TABLE 9 

Fit to Double Cluster Test, texp = 240 ks 



Type 


r-ci (kpc) 


/3i So,i 


(cts s ^cm ^arcsec 2) 


(kpc) 


02 50,2 


(cts s -'^cm 2 arcsec 2) 


xVd-o.f. 


Model 
Fitted 


50 

47 14+2-94 


1 

1.0 (F) 


3.78 X 10-8 
(3.90l°i^) X 10-8 


200 
198.00lg:7| 


0.7 
n 701+0 029 


1.38 X 10"* 
(1.40lH^) X 10-8 


213.93/193 



TABLE 10 

Spectral Fit for a Single Cluster, tgxp = 4 Ms 



Type 


Tspcc (kcV) 


^ (4;rn5^T^/"^-HrfV) (cm-5) 


xVd.o.f. 


Model 
Fitted 


4.0 
4 01+0 09 


2.547 X 10-" 
(2.55410 0^1) X 10-4 


91.40/424 



TABLE 11 

Systematic Trends in Radial Fits with Increasing 
Exposure Time 



texp (ks) 


-^single 


-'^double 


X^/d.0.f.si„gie 


X^/d.O.f. double 


240 


0/3 


0/5 


1.01 


1.11 


1000 


0/3 


1/5 


1.07 


1.20 


4000 


2/3 


5/5 


1.88 


2.20 
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Fig. 13. — Isothermal cluster temperature profile. Dashed line is input temperature 
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